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Abstract 

We examine the Mielnikov criterion for transition to chaos in case of one degree of 
freedom nonhnear oscillator with non symmetric potential. This system subjected 
to an external periodic force shows homoclinic transition from regular vibrations to 
chaos just before escape from the potential well. Especially we study the effect of 
a second resonant excitation with different phase on the system transition to chaos 
and propose a way of its control. 
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1 Introduction 



In this note shell examine vibration and show possible non feedback control 
of chaos in a simple, one degree of freedom, system subjected to external 
excitation with a non-symmetric stiffness given by the following equation: 

X + ax + 5x + = /i cos ut (1) 



where x is displacement ax is linear damping, jj, cos ujt is an external excitation 
while 5x and are linear and quadratic force terms: 

F{x) = -6x - (2) 
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Such systems, where quadratic term brakes the symmetry of a potential V{x) 
V{x)^V{-x), and = (3) 

have been a subject of studies for many years [1,2,3,4,5,6,7,8,9,10]. These in- 
vestigations were motivated by possible applications in description of physical 
systems mostly mechanical [1,3,7,8] and electrical systems [4]. They were also 
linked to possible metastable states of atoms and appear in problems within 
the elastic theory [2,11]. 

Systems which show homoclinic orbits and can be tackled analytically by per- 
turbation methods. Namely by the Mielnikov method treating [12,13] excita- 
tion and damping terms in higher order. Such a treatment has been performed 
to selected problems with both symmetric and non-symmetric nonlinear forces 
[12,13,10] to derive the necessary condition for transition to chaotic motion. 
On the other regular and chaotic regions of solutions in system parameters can 
be stabilised by using an additional weak resonant excitation [14,15,16,17]. In 
this note we shell apply this method to the nonlinear system given by Eq. 1. 
Combined with the Mielnikov approach it will predict analytically the range 
of parameters which tam the chaotic behaviour. The paper is divided into 
5 sections. After this introduction (Sec. 1) we perform Mielnikov analysis in 
Sec. 2. This discussion is followed by Sec. 3 where we include a week reso- 
nant excitation term. Its useful role in system control is shown there. The 
analytic predictions are confirmed by means of numerical simulations (Sec. 4). 
We ending up in Sec. 5 with conclusions. 



2 Melnikov analysis 



Note in this section we follow the discussion initiated by Thompson [2] where 
he derived the analytic formula for a critical amplitude of a nonlinear oscillator 
described by a similar (to our Eq. 1) equation. We decided to include this 
section as an important introductory part to our main results to be given in 
the next section. 

Thus we start our study from the second order equation of motion 1. Trans- 
forming it into two differential equations of the first order we get 

x^v (4) 
V — —5x — -|- e [—av -\- fl cos ut] . 

Looking for stable and unstable manifolds we have introduced small parameter 



2 



e to the above equations and renormalised parameters a and /i via a = ea 
and fi = e/i, respectively. 




Fig. 1. (a) The potential well of the unperturbed Hamiltonian (Eq. 4) for 6 = 1 
and 7 = 1.089, the energy level V = E corresponds to the Hamiltonian system with 
a tangent point at xq ; (b) stable and unstable manifolds for unperturbed 
(gray) and a damped and excited system (in black). The distance d between 
and can be described by Mielnikov function M(fo) (Eq. 9). Note xq indicates 
the local extremum of potential V{x) (Fig. la) which simultaneously represents a 
saddle point in the phase plane (Fig. lb). 

The corresponding unperturbed Hamiltonian have the following form: 

H'='^ + V{xl (5) 

where 

VW = |^ + ^ (6) 

is the potential of a non-symmetric well plotted in Fig. la. This function has 
the local peak at the point 

^0 = --. (7) 

7 

Existence of this point with a horizontal tangent make possible homoclinic 
bifurcations of the system i.e. potential transition from a regular to chaotic 
solution. After simple integration (Appendix A) we get homoclinic orbits (Fig. 
lb) as: 
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Note the characteristic saddle point xq is going to be reached in exactly defined 
albeit infinite time t corresponding to +00 and —00 for stable and unstable 
orbits, respectively. 

In case of perturbed orbits and the distance between them is given 
by the Mielnicov function M{to): 

+00 

M{to) = J h{x*,v*) A g{x*,v*)dt (9) 
—00 



where the corresponding differential forms h as the gradient of unperturbed 
Hamiltonian (Eq. 3) leading to equations of motion 



while g as its perturbation form of the above (Eq. 4): 



h— (Sx + ^x^^ dx + vdv, (11) 
g — {p> cos CUT — av) dx 

are defined on homoclinic manifold {x,v) = {x*,v*) (Eq. 6, Fig. 2b). From the 
above (Eqs. 8-12) the Mielnikov integral is given by 

+00 

M{to) ^ j (Jlv* cos {ujt) - (12) 



After substituting x*{t) and v*{t) by formulae given in Eq. 8 and taking r = 
\/5t/2 we get: 

^(^o) = - 3/i-/2, (13) 

z 7 7 



where 



0.0 0.4 0.8 1.2 1.6 -1.0 -0.5 0.0 0.5 1.0 

(a) 0) (b) >p/^ 

Fig. 2. (a) Critical value of fic plotted against u for system parameters: 7 = 1.089, 
6 = 1.0. Control parameters were taken as ^ = vr while F = —0.45, 0.0 and 0.45. 
The singular point below the upper curve corresponds to /i = 0.1 and uj = 0.85 
used in numerical calculations - Figs. 4-5. (b) Phase diagram in a ^-F plane for 
a number of /x values, corresponding lines. The arrow indicate the direction of 
increase: ^ = 0.04, 0.05, 0.06, 0.08, 0.10 and 0.15, respectively. Here 'R' and 'C 
denotes regions with regular and chaotic solutions. 



00 

dt 



tanh^ T 
cosh^ T 



—00 
+00 



—00 



tanhr 2uj (r + tq) 

-o— cos l=r- 

cosh r 



i2= I dt — -y- cos 7= (14) 



After evaluation of these elementary integrals we get condition of homoclinic 
transition to chaos, corresponding to a horseshoe type of cross-section and can 
be written as: 

V MW = 0. ^/O. (15) 



Evaluating above integrals (Eq. 14) after some lengthly algebra the last con- 
dition (Eq. 15) leads to a critical value of excitation amplitude /Xc for which 

^ic = T ^smh ^ . 16 

5 7ru;^7 



As a result we get the critical amplitude /ic versus frequency u;, which is 
plotted in Fig. 2(a) (the middle curve for F = 0). One should note here that 
the function has a local minimum at the point uj ~ 0.6. Note that in spite 
of using cos(ci;t) in place of sin(ti;t), in the external excitation, discussed by 
Thompson [2] in the external excitation this result appears to be the same. 
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3 Effect of a weak resonant excitation 

Let us consider an additional excitation term in starting equation (Eq. 1) 

X + ax + 5x + ^x? = fi cos (ut) + jiF cos {ut + (17) 

where F is a scaling coefficient while ^ represents a phase of a weak resonant 
excitation. Now instead of Eq. 4 one can separate unperturbed and small 
perturbation parts in the following differential equations of the ffist order: 



x = v (18) 
i) = —5x — 'jx'^ + e [—av + jl cos (cut) + flF cos {cut + \E')] . 

Now one can repeat most of calculations from the previous section assuming 
that the excitation term is composed of two parts. After simple algebra we 
get: 

cos (cut) + F cos {ut + ^) = Vl + F2 + 2F cos ^ cos {cut + a), (19) 



where 



,'l + Fcos^\ , , 

a = arctan — — (20) 

^ Fsin^ ' ^ ^ 





Fig. 3. Function B{F,^) for ^ = 0.1 (a) and 0.05 (b) respectively. 

Therefore in this case the critical amplitude fic depends on control parameters 
F and \Ef: 



1 6^V6a , 
= 7 5— smh 

5 iru^'y 



IjJTX 



1 



+ F2 + 2Fcos^ 



(21) 
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This is our principal result in this paper. In Fig. 2a we have plotted the 
results of /i^ for = vr and F = -0.45, 0.0, 0.45. Note that in case of F = 
the additional weak resonant excitation is absent but for any other cases it 
influences the system vibrations. It can drive the system away {F = 0.45), 
or into {F = —0.45) the regions with a potential chaotic solution. To show 
the separation of regions with regular solutions and potential chaotic ones we 
calculated function B{F, ^) ~ M(to) for a properly chosen to (Eq. 15): In our 
case we defined B as follows: 



a 



B = - \ sinh I — 1 Vl + i^^ + 2Fcos^. (22) 

5 7rci;^7 



In Fig. 2b we present the corresponding nodal lines for given amplitudes /i = 
0.04, 0.05, 0.06, 0.08, 0.10 and 0.15, respectively. On the other hand in Figs. 
3a and b we show the i? as a function of F and for = 0.1 and 0.05 (3a 
and b respectively). Note, the nodal lines separating regular region from the 
regions where chaotic solutions occur can be easily obtained by cross-section 
with 5 = plane. 



(a) 




0.6 -0.2 0.2 0.6 



(b) 




0.6 -0.2 0.2 0.6 



Fig. 4. Phase portraits plotted with lines together with Poincare sections denoted 
by points for F = (a) and F = 0.45 (b). /i = 1 and all other system parameter as 
in Fig. 2a. Note system parameters correspond to the point denoted with arrow in 
Fig. 2a. 



4 Numerical results 



To confirm the analytic predictions we have done numerical simulations for 
system parameters chosen as in point denoted with arrow in Fig. 2a. In Fig. 
4 we have plotted phase portraits and corresponding Poincare sections for 
F = (a), F = 0.45 (b) and \1/ = vr. As expected from our initial discussion, 
based on approximate Mielnikov approach, we find chaotic solution in Fig. 
4a. The maximal Lyapunov exponent appeared to be positive ( Ai ~ 0.05. 
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On the other hand Fig. 4b shows a regular type of motion synchronised with 
a periodic excitation force. Note that the above numerical (Fig. 4) results 
confirms the analytic investigations (Fig. 2). To show the fractal structure the 
strange attractor (in Fig. 4a) has been also magnified in Fig. 5. Note, different 
range of axes in Fig. 5a-c. The characteristic double lines structure occur with 
the scaling factor which was estimated to be 4. 
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Fig. 5. (a) Chaotic attractor (magnified from Fig. 4 and its fractal structm'e (b-c). 



5 Last remarks and conclusions 



Using the Mielnikov method we have got the analytical formula for transition 
to chaos in a one degree of freedom, system subjected to parametric excitation 
with a non-symmetric stiffness with self-excitation term. Including a weak 
resonant term we have been able to control the type of solution. Namely it 
suppressed the chaotic solution by shrinking the lower bound lines of chaos in 
a 'phase' diagram (Fig. 2b). 

In our numerical simulations of system vibration we have used the same initial 
conditions Xini = Xq + 0.1 and Vini = 0.1 representing the point in the phase 
plane relatively close saddle one {x,v) = (a;o,0). 
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Our basic chaotic solution were very close to that one obtained by Thompson 
[2] with (I7I = 1.0 and fi = 0.1089). Although we have not discussed the 
problem of escape from the potential well we changed slightly his parameters 
(7 1.089 and // — > 0.1) purposely. This has helped us to have more robust 
chaotic solution away the escape point. 

In the present analysis we have used only positively defined 7 and 6 (see Fig. 
la). Nevertheless we claim that our results are general. Note that another 
choice of 7 sign would only mirror our solution along the x axis {x — > —x). On 
the other hand changing 5 — > —\5\ would also rescale x axis and simultaneously 
shift the saddle point as well as the local minimum of the potential: 



where x' — x — \S\/j. 
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Appendix A 

Starting from unperturbed Hamiltonian (Eq. 5) we note that at the saddle 
point Xq (Eq. 7) the system velocity is zero v = and the energy has only its 
potential part 




(23) 



E^V{x 




(A.l) 



transforming Eq. 5 for chosen constant energy (Eq. Al) 



E — -— — const 
67'' 



(A.2) 



we get the following expression for velocity: 




(A.3) 
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Now one can perform integration over x 



dx 



1 



(A.4) 




1 53 

672 



5x£ 
2 



where to represents an integration constant. Finally, we get so called homo- 
clinic orbits (Fig. lb). 
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